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Dynamics of a periodically time dependent quantum sys- 
tem is reflected in the features of the eigenstates of the Flo- 
quet operator. Of the special importance are their localiza- 
tion properties quantitatively characterized by the eigenvec- 
tor entropy, the inverse participation ratio or the eigenvector 
statistics. Since these quantities depend on the choice of the 
eigenbasis, we suggest to use the overcomplete basis of coher- 
ent states, uniquely determined by the classical phase space. 
In this way we define the mean Wehrl entropy of eigenvectors 
of the Floquet operator and demonstrate that this quantity 
is useful to describe quantum chaotic systems. 
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I. INTRODUCTION 

Analysis of quantum chaotic systems is often based on 
the statistical properties of the spectrum of the Hamiho- 
nian H (in the case of autonomous systems) or the Flo- 
quet operator F (in the case of periodically perturbed 
systems). In general, quantized analogous of classically 
chaotic systems display spectral fluctuations conforming 
to the predictions of random matrices. Depending on the 
geometrical properties of the system one uses orthogonal, 
unitary or symplectic ensemble [P|] . 

Another line of research deals with eigenstates of the 
analyzed quantum system. One is interested in their 
localization properties, which can be characterized by 
the eigenvector distribution |^-|^ the entropic localiza- 
tion length ITJ or the inverse participation ratio S. All 
this quantities, however, are based on the expansion of 
an eigenstate in a given basis {&i}, which may be chosen 
arbitrarily. If one chooses (with a bad will), {bi} as the 
eigenbasis of F, all these quantities carry no information 
whatsoever. One may ask, therefore, to what extend the 
quantitative analysis based on the eigenvector statistics 
is reliable. 

Let G denote a unitary operator, such that {bi} is its 
eigenbasis. We showed [p|-pl| that the eigenvector statis- 
tics of a quantum map F conforms to the prediction of 
random matrices, if operators F and G are relatively ran- 
dom, i.e., their commutators are sufficiently large. 



In this paper we advocate an alternative method of 
solving the problems with arbitrariness of the choice of 
the expansion basis. Instead of working in a finite dis- 
crete basis, we shall use the coherent states expansion 
of the eigenstates of F. For several examples of compact 
classical phase spaces one may construct a canonical fam- 
ily of the generalized coherent states ||l^. Localization 
properties of any pure quantum state may be character- 
ized by the Wehrl entropy, equal to the average log of 
its overlap with a coherent state |l^,|l^. We propose to 
describe the structure of a given Floquet operator F by 
the mean Wehrl entropy of its eigenstates. This quantity, 
explicitly defined without any arbitrariness, is shown to 
be a useful indicator of quantum chaos. 

This paper is organized as follows. In section II we 
review the definition of the Husimi distribution, stellar 
representation, and the Wehrl entropy. For concreteness 
we work with the SU{2) vector coherent states, linked to 
the algebra of the angular momentum operator and corre- 
sponding to the classical phase space isomorphic with the 
sphere. In section III we define the mean Wehrl entropy 
of eigenstates and present analytical results obtained for 
low dimensional Hilbert spaces. Exemplary application 
of this quantity to the analysis of the quantum map de- 
scribing the model of the periodically kicked top is pro- 
vided in section IV. 



II. HUSIMI DISTRIBUTION AND STELLAR 
REPRESENTATION 

Consider a compact classical phase space J7, a classi- 
cal area preserving map Q : Q —^ n and a correspond- 
ing quantum map F acting in an A'^-dimensional Hilbert 
space Ti.N- A link between classical and quantum me- 
chanics can be established via a family of generalized 
coherent states \a). For several examples of the classical 
phase spaces there exist a canonical family of coherent 
states. It forms an overcomplete basis and allows for an 
identity resolution J^ \a){a\da = 1. Any mixed quantum 
state, described by a density matrix p can be represented 
by the generalized Husimi distribution pq] , (Q-function) 



Hp{a) := {a\p\a). 



(2.1) 



The standard normalization of the coherent states, 
{a\a) — 1, assures that L^Hp{a)da = 1. For a pure 
quantum state \^) the Husimi distribution is equal to 



the overlap with a coherent state H^{a) := K-i/jja)!^. Let 
us note that the Husimi distribution was successfully ap- 
plied to study dynamical properties of quantized chaotic 
systems [||,0. 

Consider a discrete partition of the miity into n terms; 
J27=iPi = 1- The Shannon entropy Sd = -J27=iPi^'^Pi 
characterizes uniformity of this partition. In an analo- 
gous way one defines the Werhl entropy of a quantum 
state p p3| 



Hp{a) hi[Hp{a)]da, 



(2.2) 



in which the summation is replaced by the integration 
over the classical space Jl. This quantity characterizes 
the localization properties of a quantum state in the clas- 
sical phase space. It is small for coherent states localized 
in the classical phase space ft and large for the delocalized 
states. The maximal Wehrl entropy corresponds to the 
maximally mixed state p*, proportional to the identity 
matrix, for which the Husimi distribution is uniform. 

Although the notions of the generalized coherent 
states, the Husimi distribution, and the Wehrl entropy 
are well defined for several classical compact phase 
spaces, in this work we analyze in detail only the most 
important case ft = S'^. This phase space is typical to 
physical problems involving spins, due to the algebraic 
properties of the angular momentum operator J. In this 
case one uses the family of spin coherent states |z?, f) lo- 
calized at points {'d,ip) of the sphere S^. These states, 
also called SU(2) vector coherent states, were introduced 
by Radcliffe |l^ and Arecchi et al. [19l and are an exam- 
ple of the general group theoretic construction of Perelo- 
mov (l|. 

Consider an A^ = 2j -I- 1 dimensional representation of 
the angular momentum operator J. For a reference state 
one usually takes the maximal eigenstate \j, j) of the com- 
ponent Jz- This state, pointing toward the "north pole" 
of the sphere, enjoys the minimal uncertainty. The vec- 
tor coherent state represents the reference state rotated 
by the angles 'd and ip. Its expansion in the basis |j, m-), 
TO = -j, . . . , +j reads @ 

\,9,^)= ^ sin-'-"(|)cos^+'"(|)x 



exp 



(*(J 



— m)(p 



2j 
j — m 



1/2, 



I j, to) 



(2.3) 



where /p dip J^ sinMi3|i?, (^)(t?, ip\{2j + \)/Aii = 1. 

For the SU{2) coherent state the distribution (2.1) 
equals in this case H^{'d, (p) :— \{ip\i9, (/?)p. Two different 
spin coherent states overlap unless they are directed into 
two antipodal points on the sphere. The Husimi repre- 
sentation of a spin coherent state has thus one zero (de- 
generated A — 1 times) localized at the antipodal point. 
In general, any pure quantum state can be uniquely de- 
scribed by the set of A — 1 points distributed over the 



sphere. Some of these zeros may be degenerated, just as 
in the case of a coherent state. This method of charac- 
terizing a pure quantum state is called the stellar repre- 
sentation |pi| , p2| . 

In the analyzed case of the classical phase s pace iso- 
morphic with the sphere S"^ the Wehrl entropy (2^) of a 
state p equals 



S„^- 



2j + l 
47r 



2ir 



sin Ml? / dipHpid, ip) In Hpi-d, ip) 



(2.4) 



since the measure element da is equal to (2j -t- 
1) sm'dd'ddip/4'K. Under this normalization the entropy 
of the maximally mixed state p, equals to In A. 

The Husimi distribution of an eigenstate |j, m) may be 
computed directly from the definition (2.3). Due to the 
rotational symmetry it does not depend on the azimuthal 
angle ip 



H, 



j-m) 



(^) := sin2(J-'")(^/2)cos2(J+™)(79/2)f 



2j 



(2.5) 



which simplifies the computation of the Wehrl entropy. 
Simple integration gives for the reference state | j, j) 



•Jcoh 



A- 1 

A 



2j 



2.? + 1 



(2.6) 



Due to the rotational invariance the Wehrl entropy is the 
same for any other coherent state. 
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Table 1. Wehrl entropy S\j^m) for the eigenstates of 
Jz and its mean Sj, for A — 2,3,4,5,6. Due to the 
geometrical symmetry S\j^m) — "5*1 j, -m)- 

The Wehrl entropies for other eigenstates of Jz are col- 
lected in Tab. 1 for some lower values of A. These results 
may be also obtained from the general formulae provided 
by Lee [g3| for the Wehrl entropy of the pure states in the 
stellar representation. Eigenstate |j, m) is represented by 
J -l-TTi zeros at the south pole and j — m zeros at the north 
poles. For j = 1/2 (A = 2) all the states are SU{2) co- 
herent, so their entropies are equal. For j ~ 1 (A ~ 3) 



the coherent state |1, 1) is characterized by the smallest 
entropy, while the state 1 1 , 0) by the largest (among the 
pure states). The larger N, the more place for a various 
behaviour of pure states, measured by the values of S. 
The axis of the Wehrl entropy is drawn schematically in 
Fig.l. 
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FIG. 1. Axis of Wehrl entropy for pure states of N dimen- 
sional Hilbert space; a.) N = 2 for which Smin ~ 5max = 1/2; 
b) iV = 3 for which Smin = 2/3, % ^ 0.77, {8)3 ^ 0.83, 
S'max « 0.973; c) N = 4 for which Smin = 3/4, Sj, « 0.95, 
{8)4 ~ 1.08, Smax ~ 1.24; and d) iV >> 1, where 
Smin = 1 - 1/N while {S)n ~ InA^ - 0.423. 

It has been conjectured by Licb ||23] that vector coher- 
ent states are characterized by the minimal value of the 
Wehrl entropy S'min = •S'coh, the minimum taken over all 
mixed states. For partial results in the direction to prove 
this conjecture see P5| , p3p(^ ] . It was also conjectured 
p3[ that the states with possibly regular distribution of 
all TV — 1 zeros of the Husimi function on the sphere 
are characterized by the largest possible Wehrl entropy 
among all pure states S'max- Such a distribution of zeros 
is easy to specify for (iV — 1) = 4, 6, 8, 12, 20, which cor- 
respond to the Platonian polyhedra. For A'' = 2 all pure 
states are coherent, so S'min = Smax = 1/2- For iV = 3 
the maximal Wehrl entropy Smax = 5/3 — ln2 « 0.97 is 
achieved for the state |1,0), for which the two zeros of 
Husimi function are localized at the opposite poles of the 
sphere. For N = 4 the state with three zeros located 
at the equilateral triangle inscribed in a great circle is 
characterized by Smax = 21/8 - 2 In 2 w 1.24. It wiU be 
interesting to find such maximally delocalized pure states 



for larger values of N, and to study the dependence Smax 
of TV. 

Let us emphasize that for iV >> 1 the pure states ex- 
hibiting small Wehrl entropy, (of the order of Smin), are 
not typical. In the stellar representation coherent states 
correspond to coalescence of all iV — 1 zeros of Husimi 
distribution in one point. In a typical situation the den- 
sity of the zeros is close to uniform on the sphere, and the 
Wehrl entropy of such delocalized pure states is large. A 
random state can be generated according to the natural 
uniform masure on the space of pure states by taking any 
vector oi a, N X N random matrix distributed according 
to the Haar measure on U{N). Averaging over this mea- 
sure one may compute the mean Wehrl entropy {S)n of 
the pure states belonging to the N dimensional Hilbert 
space. Such integration was performed in P,p7 28 1 in a 
slightly different context leading to 



^ 1 



(2.7) 



where VP denotes the digamma function. Note that an- 
other normalization of the coherent states used in Ref. 
[p8[ , leads to results shifted by a constant — IniV. Such 
a normalization allows one for a direct comparison be- 
tween the entropies describing the states of various N. 
In the asymptotic limit N -^ 00 the mean entropy {S)n 
behaves as In A^ + 7 — 1 ~ In A^ — 0.42278, which is close 
to the maximal possible Wehrl entropy for mixed states 
Sp^ = IniV. This result is schematically marked in Fig 
Id. 



III. MEAN WEHRL ENTROPY OF 
EIGENSTATES OF QUANTUM MAP 

Consider a quantum pure state in the A^-dimensional 
Hilbert space. Its Wehrl entropy computed in the vector 
coherent states representation may vary from 1— 1/iV, for 
a coherent state, to the number of order of In N, for the 
typical delocalized state. This difference suggests a sim- 
ple measure of localization of cigcnstates of a quantum 
map F. Denoting its eigenstates by |V'i); i = 1, . . . , A^ we 
define the mean Wehrl entropy of eigenstates 



Sf 



1 ^ 



i,,)- 



(3.1) 



This quantity may be straightforwardly computed nu- 
merically for an arbitrary quantum map F. For quan- 
tum analogues of classically chaotic systems exhibiting 
no time reversal symmetry all eigenstates are delocalized. 
In this case the mean Wehrl entropy of eigenvectors Sp 
fluctuates around {S)n ^ InN. 

In the opposite case of an integrable dynamics the 
eigenstates are, at least partially, localized. A simple 
example is provided by any Hamiltonian diagonal in the 



Jz basis (or the basis of any other component of J). The 
mean Werhl entropy Sj^ is given in table 1 for some val- 
ues of N . Further analysis shows ^1|] that for larger N 
the mean entropy behaves as ^IniV. This result has a 
simple interpretation. Let us divide the surface of the 
sphere into N ^ hT cells. A typical eigenstate of Jz is 
localized in a longitudinal strip of a constant polar angle 
9 and covers \/N of the cells, so its entropy is of the order 
of In y/N. 

The quantity S is well-defined in a generic case of op- 
erators F with a nondegenerate spectrum. In the case 
of the degeneracy, there exists a freedom of choosing the 
eigenvectors; to cure this lack of uniqueness we define Sp 
as the minimum over all possible sets of eigenvectors of 
F. Having a general definition of the mean Wehrl entropy 
of eigenvectors of an arbitrary unitary operator, one may 
pose a question, for which operators F,„i,i (F,„ax) of a 
fixed N with a nondegenerate spectrum this quantity is 
the smallest (the largest). It is clear that Sp^i^ is larger 
than S'coh (for A^ > 2), since the set of any A'' coher- 
ent states does not form an orthogonal basis. On the 
other hand, the minimum is smaller than 5j, , as explic- 
itly demonstrated in Appendix for A^ = 3. The value 
•S^i^max i^ larger than the average over the random unitary 
matrices {S)ij{n) = {S)n and smaller than Sp^ = InN. 

The mean Werhl entropy of the eigenstates Sp may 
be related with the eigenvector statistics of the oper- 
ator F. Let us expand a given coherent state in the 
eigenbasis of the Floquet operator, |a) = J2i=i Ci(ck)IV'i)- 
The dynamical properties of a quantum system are char- 
acterized locally [g9[ by the Shannon entropy Ss{a) :— 
— J2i=i |ci(a)P In |ci(a)p. The mean Wehrl entropy may 
be thus written as an average over the phase space 



Sf = — Ss{a)da. 



(3.2) 



This link is particularly useful to analyze the influence 
of the time reversal symmetry. In presence of any gen- 
eralized antiunitary symmetry the operator F may be 
described by the circular orthogonal ensemble (COE). 
There exists a symmetry line in the phase space and the 
coherent states located along this line display eigenvec- 
tor statistics typical of COE |3^] . This symmetry is also 
visible in the stellar representation of the eigenstates and 
manifests itself by a clustering of zeros of Husimi func- 
tions |^,Q. However, a typical coherent state does not 
exhibit such a symmetry and its eigenvector statistics 
is typical to the circular unitary ensemble (CUE). Thus 
for a system with the time-reversal symmetry the mean 
Wehrl entropy will be slightly smaller than for the analo- 
gous system with the time reversal symmetry broken, but 
much larger than the Shannon entropy of real eigenvec- 
tors of matrices pertaining to the orthogonal ensemble. 



IV. MEAN WEHRL ENTROPY FOR THE 
KICKED TOP 

In order to demonstrate usefulness of the mean Wehrl 
entropy in the analysis of quantum chaotic systems we 
present numerical results obtained for the periodically 
kicked top. This model is very suitable for investigation 
of quantum chaos [p5[[l|| . Classical dynamics takes place 
on the sphere, while the quantum map is defined in terms 
of the components of the angular momentum operator J. 
The size of the Hilbert space is determined by the quan- 
tum number j and equals A^ — 2j + 1. One step evo- 
lution operator reads Fq = exp(— zpj^) exp(— i/cJj/2j). 
For p = 1.7 the classical system becomes fully chaotic for 
the kicking strength fc « 6 p^ . This system possesses a 
generalized antiunitary symmetry and can be described 
by the orthogonal ensemble. The time reversal symme- 
try may be broken by an additional kick M . The system 
Fu = Foexp{—ik'Jy/2j) pertains to CUE and will be 
called the unitary top. 




FIG. 2. Husimi distribution of exemplary eigenstates of the 
Floquet operator of the orthogonal kicked top for A = 62 
in the dominantly regular regime (A: = 0.5), a) and b), and 
chaotic regime (fc — 8.0), c) and d). The sphere is represented 
in a rectangular projection with t — cos §. 

Fig. 2 presents the Husimi distributions of two eigen- 
states of Fo for p = 1.7 in the regime of regular motion 
{k — 0.5) and two, for which the classical dynamics is 
chaotic (fc — 8.0). In the quasiregular case the eigen- 
states are localized close to parallel strips, covered uni- 
formly by the eigenstates of Jz- On the other hand, the 
eigenstates of the chaotic map are delocalized at the en- 
tire sphere. These differences are well characterized by 
the values of the Wehrl entropies, equal correspondingly: 
a) 2.77, b) 2.66; and c) 3.72, d) 3.80. The data, obtained 
for A^ = 62, may be compared with the mean entropy of 



the unperturbed system, Sj^ « 2.465, the mean Wehrl 
entropy of chaotic system without time reversal symetry, 
{S)q2 ~ 3.712, and the maximal entropy of the mixed 
state, Sp, = ln62 « 4.1271. 

The above eigenstates are typical for both systems, and 
the other 60 states display a similar character. The prop- 
erties of all eigenstates are thus described by the mean 
Wehrl entropy of eigenstates Sp ■ The dependence of this 
quantity on the kicking strength k is presented in Fig. 3. 
To show a relative difference between the entropies typi- 
cal to the regular dynamics we use the scaled coefficient 



f,{F) := 



Sf — Sj^ 



N 



(4.1) 



Per definition fi is equal to zero, if F is diagonal in the 
Jz basis, which corresponds to the integrability. In the 
chaotic regime F is well described by CUE and n is close 
to unity. This is indeed the case for the unitary top 
with k' = fc/2 and A: > 6. The growth of /i is bounded 
and therefore it cannot, in general, follow the increase of 
the classical Kolmogorov-Sinai entropy A (the Lapunov 
exponent averaged over the phase space), which for the 
classical system grows with the kicking strength k pi. 
The data for the orthogonal top fluctuate below unity, 
due to existence of the symmetry line. The difference 
between the coefficients /i obtained for both models does 
not depend on the kicking strength, but decreases with 
N and vanish in the classical limit N ^ oo. 
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V. CONCLUDING REMARKS 

The Wehrl entropy of a given state characterizes its 
localization in the classical phase space. We have shown 
that the mean Wehrl entropy Sp of eigenstates of a given 
evolution operator F may serve as a useful indicator of 
quantum chaos. Let us emphasize that this quantity, 
linked to the classical phase space by a family of coher- 
ent states, does not depend on the choice of basis. This 
contrasts the others quantities, like eigenvector statis- 
tics, localization entropy, inverse participation ratio, of- 
ten used to study the properties of eigenvectors. It will 
be interesting to find the unitary operators (or rather the 
repers) for which Sp is the smallest or the largest. 

The mean Wehrl entropy of eigenstates enables one 
to detect the transition from regular motion to chaotic 
dynamics. On the other hand, it is not related to the 
classical Kolmogorov-Sinai entropy (or to the Lapunov 
exponent), so it cannot be used to measure the degree of 
chaos in quantum systems. Such a link with the classical 
dynamics is established for the coherent states dynam- 
ical entropy of a given quantum map |36| , |28|] , but this 
quantity is much more difficult to calculate. Both these 
quantities characterize the global dynamical properties of 
a quantum system, in contrast to the entropy of Mirbach 
and Korsch |^,|8|, which describes the local features. 

Mean Wehrl entropy characterizes the structure of 
eigenvectors of F, and is not related at all to the spec- 
trum of this operator. Thus it is possible to construct 
a unitary operator with a Poissonian spectrum and the 
delocalized eigenvectors. Or conversely, one may find an 
operator with spectrum typical to CUE and all eigen- 
states localized. This shows that the relevant informa- 
tion concerning the dynamical properties of a quantum 
system described by an unitary evolution operator F is 
contained as well in its spectrum and in its eigenstates. 
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sions and a constant interest in the progress of this re- 
search. I am also thankful to M. Kus and P. Pakoiiski 
for helpful remarks. It is a pleasure to thank Bernhard 
Mehlig for the invitation to Dresden and the Center for 
Complex Systems for a support during the workshop. Fi- 
nancial support from Polski Komitet Badah Naukowych 
in Warsaw under the grant no 2P-03B/00915 is gratefully 
acknowledged. 



FIG. 3. Scaled mean Wehrl entropy /i of the eigenvectors 
of the Floquet operator for the kicked top as a function of 
the kicking strength k for A'^ = 62. The data are obtained 
for two models: unitary top (•) and orthogonal top (o). The 
crosses denote values of the classical Kolmogorov-Sinai en- 
tropy A, which characterize the transition to chaos in the 
classical analogue of the orthogonal top. 



APPENDIX A: MINIMAL MEAN WERHL 
ENTROPY FOR TV = 3 

In the case j = 1 (A^ = 3) any pure state may be 
described by the position of two points on the sphere - 
zeros of the corresponding Husimi distribution. In the 
stellar representation the eigenstates \m) of Jz are de- 
scribed by 1 -I- m zeros at the south pole of the sphere, 
and 1 — m zeros at the north pole, with m ~ —1,0, 1. 
The mean Werhl entropy of eigenstates of Jz is equal to 



1 - (ln2)/3 « 0.769. We find another set of three or- 
thogonal pure states, characterized by a smaller value of 
S, by allowing to move one of the zeros of the Husinii 
representation along the meridian ip — (and ip — n). In 
other words, let us define two orthogonal states 



|X+) :=cosx|l) +sinx|0), 
\X-) ■■= -sinx|l) +cosx|0). 



(Al) 



Keeping the state | — 1) fixed we define thus a one param- 
eter family of orthogonal basis O^ = {| — 1), \x-), lx+)}- 
The Husimi distributions of the state \x-) has two zeros 
at the polar angles {6'_,7r}, while the zeros of |x+) are 
located at {n, 27r - 0+}. Here 6'_ = (2 - c^)/{2 + c^) and 
9+ = (2c2 - l)/(2c2 + 1) with c = tanx- The angle 9 be- 
longs to [0, 27r) and the 6 = tt represents the south pole, 
while 9 = (or 9 = 2tt) denotes the north pole. The 
location of zeros of Husinii representation of eigenstates 
of Jz and 0^/4 is presented in Fig. 4. 

a) M M b) 
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FIG. 4. Cross-section of the Bloch sphere along the merid- 
ian = (0 = tt). a) zeros of the eigenstates of the operator 
Jz- I — 1) (•), |0) (*), and |1) (x); b) zeros of the eigenstates 
of 0^4: 1-1) (•)> IX-> (*), and|x+> (x). 

In the paper of Lee |^ one finds the Wchrl entropy 
of a TV = 3 state expressed as a function of the angle ui 
between both zeros of the Husinii distribution: 



2/3 + 0-/6 



a/2 'K^-? 



(A2) 



where a = sin {uj/2). The equivalent result was ear- 
lier established by Scutaru H, 5 = 2/3 + [a - ln(l + a)], 
where a = cr/(2 — cr). Substituting for uj the angles ■k~'9± 
we get explicit formulae for the Wehrl entropy of each 
state |x±)5 which allow us to write the mean entropy 
So as a function of the parameter x- The minimum is 
achieved for x = 7r/4 (so c = 1 and cos 9- cos 0+ = 1/3) 
and reads Sq^,^ = 1 - [ln(9/4)]/3 w 0.730. This results 
shows that the eigenbasis of Jz does not provide the reper 
characterized by the minimal mean Wehrl entropy, but 



it does not solve the problem of finding the global mini- 
mum. One can expect, that this minimum will occur for 
a set of N mutually orthogonal states, each of them as 
close to the coherent state, as possible. 
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